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Abstract 



Flight test and theoretical aerodynamic data were obtained for a 
flight test fixture mounted on the underside of an F-104G 
aircraft. The theoretical data were generated using two codes, a 
two-dimensional transonic code called Code H, and a 
three-dimensional subsonic and supersonic code called wing-body. 
Pressure distributions generated by the codes for the flight test 
fixture as well as boundary layer displacement thickness 
generated by the two-dimensi onal code were compared to the flight 
test data. The two-dimensional code pressure distributions 
compared well except at the minimum pressure point and trailing 
edge. Shock locations compared well except at high transonic 
speeds. The three-dimensional code pressure distributions 
compared well except at the trailing edge of the flight test 
fixture. The two-dimensional code does not predict displacement 
thickness of the flight test fixture well. 



ii 





CONTENTS 

Page 

INTRODUCTION 1 

SYMBOLS ............. 2 

FACILITY DESCRIPTION .... 3 

Flight Teat Fixture 3 

ANALYSIS CODE DESCRIPTION ... ...... 5 

Code H . . . 5 

The Wing-Body Code 6 

MODELING OF THE FLIGHT TEST FIXTURE 8 

Code H . 8 

Wing-Body Code Modeling .......... 10 

DISCUSSION AND RESULTS 10 

Code H Pressure Distributions ... ....... 10 

Wing-Body Pressure Distributions ...... .. 13 

Displacement Thickness Distributions ..... . 15 

CONCLUSIONS 16 

Code H Comparisons 17 

Wing-Body Comparisons .. ...... 17 

Displacement Thickness Comparisons ............. 18 

REFERENCES 19 



INTRODUCTIOH 



The use of theoretical prediction techniques can be & useful tool 
in most engineering applications. In the case of aerodynamics, many 
computer codes exist that aid the engineer with the design and analy- 
sis of aircraft and aircraft components. Two commonly used codes are 
a two-dimensional transonic analysis code developed by Bauer, 
Garabedian, JSaeson, and Korn described in reference 1, and a three- 
dimensional subsonic and supersonic wing-body analysis code developed 
by Frank Woodward and described in reference 2. Both analysis codes 
have been used successfully in predicting parameters for specific 
shapes. For example, the two dimensional code, hereafter referred to 
as Code H, has been used successfully for the prediction of supercri- 
tical airfoil characteristics and the three-dimensional code, here- 
after referred to as the wing-body code, has been used successfully to 
predict the characteristics of various wing fuselage configurations. 

Wind tunnel tests are frequently used, as a means of obtaining 
experimental data. However, the data obtained in such tests in most 
cases mur t be corrected in order to obtain the results valid for 
flight vehicles. Such limitations as scale effects due to Reynolds 
number, size limitations for models or test specimens due to test sec- 
tion dimension, improper scaling of noise or turbulence levels in the 
wind tunnel, and unreliable data near a Mach number of 1 .0 due to 
problems such as shock reflections off of the tunnel walls must be 
considered when conducting wind tunnel tests. The Dryden Flight 
Research Facility has developed an instrumented flight test fixture 
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{ FTF) that can be attached to the underside of an F-104G aircraft 
and used as a "flying wind tunnel". The FTF is essentially a low 
aspect ratio fin incorporating a wedge-shaped airfoil. 

A need exists to (1) verify Code H in order to see if it 
will accurately predict the aerodynamic parameters for shapes 
that differ from those for which it was developed, and (2) find 
an aerodynamic code which will predict the aerodynamic parameters 
for the shape used on the flight test fixture. The purpose of 
this study was to determine if these two codes could be used to 
successfully predict aerodynamic parameters for the FTF. In 
order to make this determination, the instrumented FTF was 
attached to the underside of an F-104G aircraft and flight data 
was obtained. At Mach numbers of 0.6, 0.7, 0.8, 0.85, and 0.9, 
pressure distributions and boundary layer displacement 
thicknesses were determined from the flight test data, and were 
compared to the predicted values obtained using Code H. Pressure 
distributions at Mach numbers between 0.6 and 1.4 were made and 
compared to those predicted using the wing-body code. 

SYMBOLS 

r pressure coefficient 

P 

c local streamwise chord of wing panel, cm 

M free-stream Mach number 

p static pressure, N/M 2 

q free-stream dynamic pressure, 0.7 M P, N/m 

R Reynolds number 
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x/c ratio of distance from leading edge to local chord length 

a angle of attack , deg 

$ angle of sideslip, deg 

6* boundary layer displacement thickness, cm 

0 boundary layer momentum thickness, cm 

FACILITY DESCRIPTION 

Flight Test Fixture 

The FTF is a low-aspect ratio, fin-like shape and is mounted on 
the underside of an F-104G aircraft. It is oriented so that the 
longitudinal axis is aligned on the aircraft's lower fuselage center 
line (see fig. "< } . It is made primarily of aluminum and weighs 
approximately 136 kg (300 lbs), has a chord length of 203 ca (80 in), 
a span of 61 cm (24 in), and a constant thickness of 1 6„3 ca (6.4 in) 
except for the forebody (see fig. 2). Two options are available for 
forebody shapes: (1) the basic FTF shcipa with a sharp leading edge 

(wedge forebody), and (2) the radiused forebocy incorporating the 
front portion of a symmetrical supercritical airfoil. Only the wedged 
forebody wa3 used in this study. The fin has its own air data system 
v;hich consists of a pitot static probe that is mounted on a boom 
extending forward from the FTF. The probe is used to measure Mach 
number, altitude, and dynamic pressure. The FTF is equipped with flush 
static pressure orifices for measurements of chorawise and spanwise 
pressure distributions, and boundary layer rakes for measurement of 





the boundary layer velocity profile. For this study, 20 static orifi- 
ces were used located on both sides of the FTF as shown in figure 3. 
Sixteen of the orifices were placed along the chord at approximately 
the 50$ span position. Four orifices were located along a spanwise 
direction shown in figure 3 for the purpose of determining spanwise 
flow conditions o Hie boundary layer rakes were mounted on both sides 
of the FTP at the approximately 90% chord and 50% semispan positions. 

A pulse code modulation (PCM) system is used for data acquisition. 
Data froa the PCM, which is capable of multiplexing 40 channels tit a 
maximum frequency of 80 Hz, is both transmitted to the ground via 
telemetry and recorded on board. All pressure measurements were 
obtained by a 48-port scan! valve and two other individual differential 
pressure transducers. ®ie pressures measured by the scan! valve and 
transducers ware referenced to the FTF boca-static pressure. 

The F-104G aircraft has its own independent instrument system and 
an aircraft flight trajectory guidance system. The trajectory 
guidance system uplinks engineering parameters calculated on a ground 
based computer to a cockpit display in real time. From this display, 
bank error during constant Mach, a and altitude turns, Reynolds 
number error, sideslip error, and Mach number error can be determined 
in real time by the pilot. A more complete description of the FTF is 


available in reference 3 
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ANALYSIS CODE DESCRIPTION 
Code H 

Francis Bauer, Paul Garabedian, David Korn, and Anthony 
Jameson from the Courant Institute of Mathematical Sciences of 
New York University, have developed a technique of computing 
supercritical airfoil sections and determining the off-design 
flow conditions. The equations of motion used in this method are 
the equations of potential flow. The flow is assumed to be 
transonic, steady, i rrotational , inviscid, compressible, and 
two-dimensional. Instead of solving the problem of computing 
shock free transonic flow past a given profile, the inverse 
problem was solved. That is, they assumed smooth transonic flow 
and then found the body which generated it. This approach was 
taken in order to eliminate certain mathematical difficulties. 

The problem was formulated by writing the equations of motion of 
the inverse problem in matrix form, extending all of the 
variables into the complex domain, and introducing characteristic- 
coordinates and then expressing the equations of motion in 
characteristic form. A treatment of compressibility is made by 
combining a regular solution with a singular solution that is 
related to the fundamental solution in the hodograph plane. The 
formulated equations are solved numerically using a finite 
difference scheme. Not only was the 
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inverse problem of determining a shape that would result in the smooth 
transonic flow treated, the off-design conditions (at different angles 
of attack and free-straam velocity) were also solved. The result was 
several computer programs which were designated Programs A through E. 

A description of the theory and programs is available in reference 1 . 

The authors of these programs have modified and improved their 
original work by introducing a better model of the trailing edge, and 
using a rotated finite different scheme that enables them to use an 
arbitrary cur va linear coordinate system. The use of an. arbitrary cur- 
valinear coordinate system permits the handling of supersonic and sub- 
sonic free -stream Mach numbers and the capturing of shock waves as far 
back on the airfoil as desired. The turbulent boundary layer is 
treated using a sesi -empirical Eethcsd and the effects of displacement 
thickness on airfoil shape is accounted for. Shock waves are handled 
by calculating weak solutions to the applicable partial differential 
equations that include one or more shock waves that satisfy an entropy 
inequality. These modifications are included in a new program 
designated Code H. A description of this program can be found in 
reference 1 . The author claims these programs provide a physically 
adequate computer simulation of the compressible potential problems of 
transonic flow for a smooth 2D shape. 

The Wing-Body Code 

The Boeing Company Uiider contract with NASA/Ames Research Center 
has developed a three-dimensional wing-body constant pressure panel 
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code for subsonic an d supersonic potential flows. The program calcu- 
lates steady pressure distributions on wing and wing-body combinations 
of arbitrary plan fora in subsonic and supersonic flow. The surface 
pressures are integrated to give the lift, drag, and pitching moment. 
The yawing and rolling moments and the side force can be determined 
for unsyresaatrieal configurations; however, for this study, only the 
pressure coefficients that were predicted for the surface of the 
flight test fixture were used. In addition, the original version of 
the Ames wing-body code was modified by personnel at Ames and an 
updated version of this program was made available to NASA's Dryden 
Flight Research Facility for use in this study. 

The method divides the wing-body combination or wing alone into 
numerous constant pressure panels. A constant source distribution is 
used in the body panels and a vortex distribution is used in the wing 
and tail panels. The code arrives at analytical expressions for the 
perterbation velocities induced at each panel. Then, the pressure 
coefficient at the panel control points are calculated in terms of the 
perterbation velocities. The forces and moments acting on the wing- 
body combination can be calculated by using a numerical integrating 
scheme. 

A further description as well as previous utilization of the code 
at the Dryden Flight Research Facility is given in reference 4. 
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MODELING OP THE FLIGHT TEST FIXTURE 


Code H Modeling 

In all cases presented, the code was run using 0° angle 
of attack^ and four smoothing iterations of the FTF coordinates were 
made before the aerodynamic shape was conforaly mapped into the unit 
circle. The circle was overlayed with both a crude grid of 80 x 15 
mesh intervals and a finer grid of 160 x 30 mesh intervals in the 
angular and radial directions. Plow calculations and boundary layer 
corrections were computed for a maximum of 20 cycles on the crude 
grid, and a saxisua of 10 cycles on the finer grid. The convergence 
tolerance, a tolerance of the maximum velocity potential and the maxi- 
mum circulation corrections, was set at 5 x 10~®. ®ie program was ran 
until the convergence tolerance was achieved. The boundary layer 
correction option of the code was used, and the transition was set at 
the 7.5% chord position. To utilize this option, a Reynolds number 
must also be specified. In this case, Reynolds numbers of 2 x 10^ and 
14 x 10 ^ were used. 

The FTF was first modeled using 46 upper and 46 lower surface 
points with a high density of points at the wedge corner of the FTF 
located at the approximate 17% chord postion. Because of the discon- 
tinuity at the corners, the code would not complete the run. Next, a 
model with 46 upper and 46 lower surface points was attempted but the 
sharp corners of the wedge were radiused and the coordinates were run 



through a separate smoothing program before being entered into 
the two-dimensional code. Figure 4 shows a comparison at this 
model shape to the actual shape of the FTF . A comparison of the 
predicted and experimental pressure distribution is shown in 
figure 5 for the 0.7 Mach case. This figure shows that the code 
predicted pressure distribution appears to predict pressure 
coefficient levels well; however, there is a difference of 
approximately 7% chord in the chordwise location of the peak 
values of pressure coefficients. In addition, the experimental 
data indicates a trend of decreasing pressure at the trailing 
edge of the FTF that is not predicted by the code. The variation 
of the peak pressure coefficient positions between the 
experimental and the computer predicted data were believed to be 
caused by the variation of shape of the actual and computer 
models . 

In order to test this hypothesis, a model consisting of 16 
upper and 16 lower surface points was run. The model having 
fewer points allowed the smoothing subroutine internal to the 
code to have a greater effect, as shown in figure 6. The shape 
of the model is very similar to the actual FTF. Figure 7 shows 
good correspondence for the positions of the minimum pressure 
coefficients determined from experimental and theoretical data. 
The noted lack of correspondence of pressure coef f f i ci ents at the 
trailing edge was not believed to be caused by modeling and will 
be discussed in the RES n LTS AND DISCUSSION section of this 
report. 
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Wing-Body Code Modeling 

Since this code is three dimensional, it is necessary to not 
only model the FTF but also the F-104G aircraft. The wing-body 
code allows for a maximum of 100 wing and 1 °0 body panels. These 
panels were divided between the fuselage of the F-104G and the 
wing and FTF. The fuselage was modeled with 9b panels and the 
wing and the FTF used 12 panels. Thirty-two panels were used for 
the F-104G wings, and 50 ->anels were used for the FTF. Figure 8 
shows the panel breakdown of the F-104G and the FTF. The F-104G 
was modeled as a cylinder with a radius of 80 cm and a length of 
1160 cm, and a conical nose of 470 cm in length. The wings of 
the F-104G were modeled as a bi -convex surface with a thickness 
ratio of 3.36%, a semispan length of 230 cm, a sweep back of the 
quarter chord of 18. 6°, and an anhedral of 10° • ^he vertical and 
horizontal stabilizers were not included in the model in order to 
allow a greater number of panels for the FTF. 

DISCUSSION AND RESULTS 

Code H Pressure Distributions 

Figures 9(a), 9(b), 9(c), and 9(d) present pressure 
distributions for the FTF using experimentally and theoretically 
determined data for Mach numbers of 0.7, 0.8, 0.85, and 0.9. The 
experimental data is based or. Mach numbers measured with the FTF 
noseboom since only the FTF was modeled in this case. Code H was used 
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to calculate the theoretically determined pressure coefficients. 
As shown in figure 9(a) where the Mach number is 0.7, the flow 
accelerates and the pressure coefficient drops from the nose to 
about the point on the surface where the discontinuity occurs 
(approximately 17% chord). Beyond this point the flow slows as 
it changes direction and the pressure coefficient increases. 
There is good correspondence between the experimental and 
theoretical data with good correspondence between the 
experimental and theoretical data with the exception of the peak 
minimum values of pressure coefficients. Code H predicts higher 
minimum pressure coefficients than are obtained from the flight 
test data. No shock exists for this Mach number, since 
supersonic flow velocities do not occur on the FTF. 

Figures 9(b) through 9(d) show data for Mach numbers of 0.8, 
0.85, and 0.9. For these Mach numbers the flow accelerates from 
the nose to the surface discontinuity, and reaches sonic 
condition. As the flow turns through an angle of approximately 
13.05° (one-half the wedge angle) at t»e discontinuity point, it 
accelerates, reaching a peak value and goes through a normal 
shock. The shock causes a rapid increase in pressure coefficient 
and a slowing of the flow to sonic velocities. These figures 
reveal that a fairly good correspondence between experimental and 
theoretical data exists; however, the peak minimum values of 
pressure coefficients differ, and as the Mach number increases 
from 0.8 to 0.85 the predicted shock location tends to shift 
beyond the 20% chord position. At 0.9. Mach number, the shock 
has shifted to about the 50% chord position in the theoretical 


data • 


The shock remains near the 20% chord position for all Mach 
numbers for the experimental data. We can conclude that Code H 
does not accurately predict shock location for the wedge shaped 
FTF near Mach 1 or high transonic speeds. 

At approximately the 70% chord position on the FTF, the two 
curves deviate and the experimental data reveal a decreasing 
pressure that the theoretical method does not predict. A 
possible explanation of why the two curves deviate near the 
trailing edge follows; the FTF can be considered to be an 
aft-facing step of height equal to one-half the width of the FTF, 
and flight measured pressure characteristics of aft-facing steps 
presented in reference 5 indicate that the base pressure does 
indeed have an effect on the pressure measured upstream of the 
aft-facing step. But the code does not account for this because 
a trailing edge of finite thickness is advanced linearly until is 
closes or exceeds chord length, whichever occurs first. 

Therefore, an attempt was made to alter the trailing edge of the 
FTF model. The model was geometrically scaled down and a 
boattai 1 was added in order to effectively accelerate the flow at 
points near the trailing edge. 

The pressure distribution for the boattai 1 FTF in figure 10 
shows that decreasing pressure at the tailing edge is evident but 
not nearly of sufficient quantity to match the experimental data. 
The boattai 1 of greater curvature was attempted, but the code 
would not run because points spaced too closely together at the 
trailing edge led to computational difficulties. Even with fewer 
points at the trailing edge, a run could not be completed because 

of the amount of curvature needed to simulate the flow at the 
trailing edge. 
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The Wing-Body Pressure Distributions 

The wing-body code was run for the F-104G/FTF model at 0.6, 
0.7, 0.8, 0.9, 1.2, and 1.4 Mach numbers. Angles of attack of 
0°, 2° and 4° were used for the 0.6 and 1.4 Mach number cases. 
These two cases which are the extremes of the Mach number range 
tested were run at different angles of attack to determine the 
effects of aircraft angle of attack on the data. Figures 11(a) 
and 11(f) show data for these cases with the wing-body code data 
for the three angles of attack and the flight test superimposed. 
Little difference exjsts for the three sets of angle of attack 
data. The data varies little wi th angle of attack. Therefore, 
the Mach number cases of 0 . 7 , 0 . 8 , 0.9 and 1.2 were run at 2° 
angle of attack which approximates the F-104's angle of attack 
during the test flights. All cases were run at 0° sideslip angle 
and the flight test data wore recorded at very small sideslip 
angles. The ptossun: coefficients given in the output of the 
wi ng-bodv code act at the centroid of the panel and represent the 
average pressure over the panel. Since the wing-body code is 
three dimensional and the model includes the aircraft, the flight, 
test data were based on at Mach numbers measured by the aircraft 
instrumentation system instead of the flight test fixture's air 
data system. 

Figures 11(a) and 11(b) depict data for Mach numbers of 0.6 
and 0.7. These figures reveal an accelerating flow from the nose 
of the flight test fixture to the surface discontinuity and a 
slowing of the flow aft of the forebody with an accompanying 
increase of pressure coefficient. There is a fairly good match 
of the flight test and 
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theoretically predicted data. The peak negative pressure 
coefficient for the flight test data is more negative than the 
predicted value, the predicted values of pressure coefficient are 
generally higher than values determined from flight test data, 
and the trend of decreasing pressure coefficients exists near the 
trailing edge of the flight test fixture for the flight test data 
that is not present for the theoretically predicted data. 

Figures 11(c) and 11(d) show c'ata for the other two subsonic 
cases of 0.8 and 0,9 Mach numbers. As discussed earlier, at 
these flight speeds the flow is accelerated by the forebody to 
supersonic speeds and a normal shock wave forms. These shock 
waves are evident in the flight test data for Mach numbers of 0.8 
and 0.9. While the shock waves are not predicted by the 
wing-body code, the point of minimum pressure coefficient does 
occur at the 20% chord position for both sets of data. This 
method is no capable of transonic shock wave prediction. In 
addition, the data indicate an inability of the wing-body code to 
predict the decreasing values of pressure coefficient near the 
trailing edge of the FTF for similar reasons as Code H results. 
However, for the 0.9 Mach number case the two sets of data 
correspond very well over a large range of chord positions. 

The data for the supersonic Mach number cases of 1.2 and 1.4 
Mach number are shown in figures 11(e) and 11(f). For both cases 
the flow over the FTF is subsonic since the one-half wedge angle 
of approximately 13.05° is large enough to cause the shock to 
detach from the nose. 
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The portion of the shock forward of the nose is normal and the normal 
shock wave creates subsonic flew over the nose of the FTF and large 
positive pressure coefficients « The flow accelerates over the fore- 
body, does not reach sonic condi tions, and decelerates from the sur- 
face discontinuity to the trailing edge. The two sets of data match 
very well over the entire FTF; however, in the 1.2 Mach number case, 
the theoretical data reflects very sharp changes of pressure coef- 
ficient about the 50^ chord position. It is possible that the wing- 
body code is introducing fuselage effects that do not occur in the 
experimental data at this particular Hach number. To check this 
possibility, the FTF was modeled excluding the F-104G aircraft. 

Figure 12 shows that the pressure distribution smooths out and com- 
pares well to the experimental data, suggesting that the code was 
indeed introducing inaccurate fuselage effects. 

It is interesting to note that the trailing edge divergence noted 
in the subsonic data changes for the supersonic data. The experimen- 
tal data indicate an increasing and then decreasing trend of pressure 
coefficients near the trailing edge of the FTF. 

Displacement Thickness Distribution 

Code K provides a semi -empirical turbulent bcur.cary-iayer correc- 
tion in the transonic flow analysis. Displacement thickness (6*) is 
calculated by relating momentum thickness and shape factor where 
momentum thickness is determined using Von Karman ‘ S equation and the 


15 




shape factor is determined semi-empirically . Because the laminar 
portion of the boundary layer is considerably smaller than the 
turbulent portion, it is not considered in the boundary-layer 
correction calculations. For the boundary-layer correction, a 
transition point must be specified. A transition location of 
7.5% chord was used in all cases. This most closely approximated 
where transition v;as thought to occur. 

In figure 13, experimental and theoretical displacement 
thickness at 85% chord are plotted vs Mach number. Reynolds 
numbers of 2 x 10 6 and 14 x 10® were used. The shapes of the 
curve agree but they are displaced. Also, the experimental data 
with approximate Reynolds numbers of 20 x 10® was quite scattered 
around 0.8 Mach. This scatter in the experimental data is 
probably caused in part by separated flow due to a shock wave 
that is reattaching. We can conclude from the data shown in 
figure 13 that the semi-empirical boundary layer used in Code H 
does not permit precise determination of the displacement 
thickness for the FTF with the wedge shaped nose. 

CONCLUSIONS 

An F-104G aircraft with an attached FTF with a wedged shaped 
forebody has been tested at NASA/Dryder. Flight Research Facility. 
Pressure distributions and displacement thicknesses have been 
determined from the flight test data. Two theoretical prediction 
methods have been used to predict similar data for the FTF at the 
flight test speeds. One is a two-dimensional method and has been 
designated Code H by Bauer, et al, who are the authors of this 
method. The other 
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ia a three-diEsnsionai method and has been designated as the wing-body 
code by its author, Frank Woodward * The results obtained from the 
comparison of flight test data with the data predicted by the codes 
follows. 

Code H Comparisons 

(1) For subsonic flight speeds and flow over the FTP, Code E ade- 
quately predicts values of pressure coefficients except at the minimum 
pressure point and at the trailing edge. Code H predicts higher 
values at the minimum pressure point and the experimental data reveal 
an increasing pressure coefficient that Code H does not predict. 

(2) For subsonic flight speeds and supersonic flow at some point 
on the FTP, the shock wave that forms is located at the approximate 
20^ chord position. Code E predicts a shifting position of the chock 
waves with increasing speed and does not adequately predict the 
increasing pressure coefficient divergence at the trailing edge of the 
FTF. Otherwise Code K adequately predicts the level of the pressure 
coefficients at other positions on the FTF. 

V?ing-Dody Code Comparisons 

(1) For subsonic speeds and flow over the FTF, the win T-body code 
adequately predicts levels of pressure coefficients except at the 
trailing edge. The decreasing pressure coefficient divergence at the 
trailing edge is not predicted by the wing-body coda. 
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(2) For subsonic flight speeds and supersonic flow at sob® point 
on the FTF, the wing-body cods which is incapable of shock wave pre- 
diction adequately predicts pressure coefficient levels except at the 
trailing edge and at the taininum pressure point. 

(3) For supersonic flight speeds and subsonic flow over the FTF 
the wing-body code adequately predicts levels of pressure 
coefficients. 

Displacement Thickness Comparisons 

The seal -empirical boundary layer used in Code H does not precisely 
predict the displacement thickness of the FTF for the two Reynolds 
numbers tasted. 
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Fig. 6 Comparison of 32 point model 
to actual shape of the FTF. 



Fig. 7 Pressure distribution 
oene rated by code H for 32 point 
model of FTF. M=0.S„ 
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Fig. 8 Panel break-down of F-104/FTF 
model used for a wing-body code. 



Fig. 9 Pressure distribution 
generated using code H. 


i 






PrssEiire cos! (Iclani (Cp) Prsssif’a iCM.'flc'tnl (Cp) 


tft-fiSflflt 

O K=.SS3 o=4.1S 
0= - .14 



bvfiigfS 



(a) H=0 .6 . 


(b) f^O.7. 


toC£f$A 

O 64= .SSiS <r«&2 


|J= — JSJ7 



Q. 

O 


5 

6 


-1.4 

- 1.2 

- 1.0 

-.8 

-.S 

-.4 

-Jt 

0 

.2 

.4 

.6 

.8 

1.0 


r 

r 


ln-Wigh« 

O MasJKK c.-= 2J2 

Pm JQ? 
Wlrsg.&.-.Sy 

O 0=2” 
<5=0" 



Shock 


tSXrQ^ - ©' 


J3SsSS^-8— O- 


§> 


.2 .4 .e .8 

Percar.t chord (x/c) 


(c) H=0 .8 . 


(d) H=0 .S . 


Fig. 11 Pressure distribution 
generated by wing-body code. 








Fig. 13 Di spl acament thickness generated 
by code H as a function of ^ach number. 
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